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Abstract 

The motion of an incompressible fluid in Lagrangian coordinates involves 
infinitely many symmetries generated by the left Lie algebra of group of 
volume preserving diffeomorphisms of the three dimensional domain occupied 
by the fluid. Utilizing a 1 + 3-dimensional Hamiltonian setting an explicit 
realization of this symmetry algebra and the related Lagrangian and Eulerian 
conservation laws are constructed recursively. Their Lie algebraic structures 
are inherited from the same construction. The laws of general vorticity and 
helicity conservations are formulated globally in terms of invariant differential 
forms of the velocity field. 
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1 Introduction 



The configuration space of an incompressible fluid is the group of volume 
preserving transformations of the three dimensional region in R 3 containing 
the fluid. The motion is generated by the left action of the group by com- 
position and hence the velocity field is right invariant. Thus, the generators 
of the right action which form the infinite dimensional left Lie algebra are 
infinitesimal symmetries of the velocity field In the fluid mechanical 

context, this is known as the particle relabelling symmetry P],Pj, also named 
as gauge transformations in Ref.|| and as trivial displacements in Ref.|J. 
It has been the subject of many investigations to discover the corresponding 
Lagrangian and Eulerian invariants in their most general form and to express 
for them the so-called general vorticity conservation law. (see Ref. [§J for a 
discussion and references). 

One the best illustrative example for this general conservation principle 
may be the equivalences of the preservation of coadjoint orbits of vorticity, 
the right invariance on the cotangent bundle of the group of volume preserv- 
ing diffeomorphisms, the Helmholtz' vorticity and the Kelvin's circulation 
theorems (see |J or section ([675| )). Another formulation of the same principle 
results, in addition, to the Eulerian conservation laws for helicity M. Both 
examples have in common the property of being derived from the particle 
relabelling symmetry. 

In spite of the fact that the particle relabelling symmetry is a direct 
consequence of the description of motion itself, many investigations on the 
structure of symmetries and invariants employ the techniques of analysing 
the defining equations for them thereby obtaining an algebraically and quan- 
titatively incomplete picture of symmetries, invariants and their connections. 

In this work we shall present a geometric framework for a systematic study 
of symmetries and invariants which will provide, by addressing explicitly to 
Lagrangian and Eulerian formulations, a better understanding of them in the 
description of motion. 

1.1 Summary and content 

We shall consider the Lagrangian description of fluid motion in the frame- 
work of a formal geometric structure to obtain a realization of the symmetry 
algebra, the corresponding Lagrangian invariants and the general vorticity 
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conservation law. 

We shall start from the observation that the Euler equations themselves 
are the conditions for the vorticity field to be an infinitesimal symmetry of the 
velocity field. We shall use this symmetry to construct a symplectic structure 
on the time-extended domain in R x R 3 . This will define the Hamiltonian 
structure of the suspended velocity field. Identifying a Hamiltonian vector 
field as a new symmetry of the velocity field, we shall be able to generate the 
infinitesimal symmetries recursively. 

The vorticity is, by construction, an automorphism of the Hamiltonian 
structure. Since any Hamiltonian vector field is also an automorphism, we 
shall conclude that the resulting infinite hierarchy of vector fields do form 
the Lie algebra of symmetries. 

Associated to each infinitesimal symmetry we shall then introduce invari- 
ant differential forms of the velocity field. Among a plethora of differential 
invariants of various degrees we shall identify a hierarchy and use it as the 
basic ingredient for the global expressions of the laws of general vorticity and 
helicity conservations. 

After a brief description of fluid motion in the next section we shall re- 
call from Ref.0 the symplectic structure of the suspended velocity field on 
the time-extended space R x R 3 . In section(0), we shall obtain the time- 
dependent symmetries on R x R 3 of the velocity field and show that an 
equivalent representation of them are Hamiltonian vector fields on R 3 . In 
section we shall give explicit expressions for the laws of general vortic- 
ity and helicity conservations. All of these will arise in connection with the 
infinitesimal automorphisms, restriction to spatial domain in R 3 , and the 
associated invariant forms of the symplectic structure. 

2 Three dimensional fluid motion 

Let the open set D G R 3 be the domain occupied initially by an incompress- 
ible fluid and x(t = 0) = xq G D be the initial position, i.e., a Lagrangian la- 
bel. For a fixed initial position xq, the Eulerian coordinates x(t) = gt{xo) de- 
fine a smooth curve in R 3 describing the evolution of fluid particles. For each 
time t G R, the volume preserving embedding g t : D — > Qt{D) = D t C R 3 
describes a configuration of fluid. A flow is then a curve t \— > g t in the space 
of all such transformations. The time-dependent Eulerian (spatial) velocity 
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field v t that generates g t is defined by 



% = ^ = (v t og t )(z )=v(t,x) (1) 

where v t o g t is the corresponding Lagrangian (material) velocity field @,[||]. 
Since g% is volume preserving, Vt(x) is a divergence-free vector field over R 3 
and Eq.([l|) is a non-autonomous dynamical system associated with it. The 
system (0) can, equivalently, be represented as an autonomous system defined 
by the suspended velocity field 

d t + v(t,x), u = v-V (2) 

on the time-extended space R x R 3 . 

The Lagrangian description of fluid motion is the description by trajecto- 
ries, that is, by solutions of non-autonomous ordinary differential equations 
flU). Equivalence of (|l|) to the autonomous system associated with (JJ) means 
that the trajectories can be obtained by describing streamlines at each time. 
In Ref . [0] , using the Eulerian dynamical equations, we constructed a formal 
symplectic structure for (g) on a time-extended domain in R x R 3 . We shall 
now summarize this construction. 



2.1 Symplectic structure 

A symplectic structure on a manifold N of even dimension 2n is 

defined by a closed, non-degenerate two-form Q. It is exact if there exists a 
one-form 6 such that Q = —d6. Darboux's theorem quarantees the existence 
of local coordinates (q % -,Pi) % = 1, ...,n in which Q has the canonical form 
dq l A dpi. The In— form (— l) n Q n /n\ is called the Liouville volume. A vector 
field V on iV is called Hamiltonian if there exists a function h on N such 
that 

i(V)(Sl)=dh (3) 

where i(V)(-) denotes the inner product with X. The identity i(V)(dh) = 
which follows from (|3]) is the expression for conservation of h under the flow 
of V. With the correspondence @ between functions and vector fields, the 
Poisson bracket of functions on iV defined by 

{f,g} = n(V f ,V a ) = n- 1 (df,dg) (4) 
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satisfies the conditions of bilinearity, skew-symmetry, the Jacobi identity and 
the Leibniz rule. This enables us to write the dynamical system associated 
with the vector field Vh in the form of Hamilton's equations 

£ = {«,*}. (5) 

Proposition 1 Let the dynamics of the velocity field v be governed by 

<9v 

— + v • Vv = F (6) 

at 

and assume that the divergence- free field B and the function if satisfy 

^-Vx(vxB) = 0, ^ + V . V ^ = (7) 

which are the frozen field equations. Then 

(1) dt + v is a Hamiltonian vector field with the symplectic two-form 

Q = -(Vv? + v x B) • dx A dt + B • (dx A dx) (8) 

and the Hamiltonian function if. 

(2) p v = — B ■ V<£> is the invariant Liouville volume density. 

(3) If moreover B = V x A for some vector potential A, then Q is exact 

Q = —d6 , -6 = ip dt + A ■ dx (9) 

where ip is determined by the equation 

dA 

— -vx(Vx A) = V(^ + ^)- (10) 

(4) In Darboux coordinates, Q = dq A dp + dt A dh can and the potential 
field has the Clebsch representation 

A = Vs + pVq (11) 

where s is the generating function of the canonical transformation. 

The Euler flow of an ideal isentropic fluid is characterized by F = — V(P/ p), 
where P and p are pressure and density, respectively. In this case, A is re- 
placed by the velocity field itself, B becomes the vorticity w = V x v and 
Eqs.(|l0D reduces to the expression 

-V = <p + P/p + v 2 /2 (12) 

for the scalar potential. The invariant volume density p^ = — w • V<p is also 
known as potential vorticity. 
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3 Symmetries 



A time-dependent vector field U = £dt + u is an infinitesimal geometric 
symmetry of dynamical system described by v if the criterion 

[d t + v, £d t + u] = (i t + v(0)(dt + v) (13) 

is satisfied. These are the most general symmetries of the system (|l|) of first 
order ordinary differential equations ||. 

The frozen field equation (|7]) for B is an expression for it to be a sym- 
metry of v. Since p^ is a conserved function, p^ x B is also a symmetry. We 
observe that this initial symmetry leaves the symplectic two-form (||) invari- 
ant. More precisely, one checks that the Lie derivative of Q with respect to 
the normalized field p~*B ■ V vanishes. 

Our intention now is to find another infinitesimal invariance of Q which is 
a symmetry of the velocity field as well. One of the best candidates for this 
is a Hamiltonian vector field because the symplectic two-form is invariant 
under the flows of Hamiltonian vector fields. This can easily be seen from 
the identity Cu = i(U)d + di{U) for the Lie derivative together with Eq.@ 
and the closure of Q. If U is such a field, called an automorphism of Q, then 
it follows from the identity 

C [p^B,U] = C p- 1 B C U - £u£ p -l B ( 14 ) 

that [p^B, U] also leaves Q invariant. Replacing p~^B with [p~ x B, U] in 
the identity ([HP we see that one can generate an infinite dimensional alge- 
bra of Hamiltonian vector fields (over simply connected domains of fluid) as 
invariants of the symplectic two-form. 

Thus, if we can find U which is also a symmetry of the velocity field, 
this so-called algebra of symplectic automorphisms of Q will be carried over 



to the symmetry algebra of the velocity field. The identity (|Tj) will then 
become the Jacobi identity of the algebra of vector fields and will enable us 
to obtain the generators recursively. 

Proposition 2 The Hamiltonian vector field 

U = p-^-B^dt + v ) + ^B + V^xVh-V] (15) 
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associated with the symplectic two-form Vt and the arbitrary smooth function 
h is an infinitesimal symmetry of v if h is invariant under the flow of v, that 
is, dh/dt = h 7 t + v(h) = 0. In this case, the brackets 

[...[[p- l B,U],U],..] (16) 

generate an infinite hierarchy of time- dependent infinitesimal Hamiltonian 
symmetries of the velocity field v. 

The infinitesimal symmetries ([IB]) are of the form of 

U k = Udt + v)+u k) fc = 0,l,2,... (17) 

where ^'s are conserved functions of v and u k s are vector fields on the 
spatial domain in R 3 . We compute £ = 0, £i = —B(h), £2 = —B 2 (h), £ 3 = 
u 2 (£i) - £ 4 = w 3 (£i) - wi(w 2 (£i)) + ul(£ 2 ), ••• to list a few of these 

functions. The vector fields u k have a well-known interpretation in symmetry 
analysis of differential equations ||. They are the unique characteristic (or 
evolutionary) forms of U k s along the velocity field v. It follows from Eq . (|i~3|) 
by direct computation that the symmetry condition for u k reduces to 

[d t + v,u k ] = 0. (18) 

We shall now switch to this equivalent representation of symmetries on 
the flow space R 3 because this form is more suitable to present their geomet- 
ric features. Since the bracket of characteristic vector fields is the same as 
the characteristic form of the bracket, the hierarchy ( fTED can be written as 
[...[[ / o~ 1 S, u], u], ...] and generates time-dependent vector fields on R 3 satis- 
fying Eqs.flnf). The Hamiltonian vector fields (0) are divergence-free with 
respect to the Liouville volume. To this end, it would be appropriate to set 
p 9 — 1. This greatly simplifies the general expressions for symmetries and 
allows a better manifestation of their generic properties. 

Using notations of three dimensional vector calculus and the divergence- 
free property of B, u 1; the hierarchy of symmetries ( |16"D in their characteristic 
form can be obtained recursively from 

u fc = (Vxii 1 xfB = Vx(u 1 xu M ), k = 2,3,4, ... (19) 

where the operator V x uix is the "curl of cross-product with ui". We find 

u x = u = V(f x Vh , u 2 = x VB(h) , 

u 3 = Vy? x V[VB(h) ■ (V<p x Vh)] , . . . , u k = V^> x Vu k ^(h) , . .(20) 
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for the first few and the generic members of the hierarchy. They are made up 
of gradients of conserved functions and this property is preserved under the 
action of symmetries. Hence, their flow lines are described by intersections 
of surfaces defined by conserved functions ip, h, B • Vy2, B • V/t, (ui • V)(B ■ 
V/i), ... of the velocity field. 

The vector fields k = 1,2,3,... are manifestly divergence-free (cf. 
Eqs.([l9l)) and are in Clebsch form (cf. Eqs.(pOl)). Then, as pointed out 
in (Exercise 1.4-1 of) Ref. 0, we can introduce vector potentials and, 
moreover, write u^'s as Hamiltonian vector fields on R 3 . From Eqs.(|i~9l) we 
can read off the potentials as 

Ai = <pVh, A k = iii x u fe _i = (u x x Vx) fe-1 A , k = 2,3,4,... (21) 

where the operator generating from the potential A of B is to "take curl 
then multiply with ui" . 

The Hamiltonian structure of m^'s is defined by the Poisson bracket 

{/ )S }, = V r V/xV 9 (22) 

which is characterized by the Hamiltonian function <p of the suspension (pj). 
The Hamiltonian functions for the vector fields ( PU| ) are 

h = h, h 2 = u (h), h 3 = -iti(u {h)) , ■■■ (23) 

from which it follows by induction that for the symmetry 

h k = (-iM-^Mh)) , k = 3,4,5, ... (24) 

is the Hamiltonian function. These are all conserved functions of the velocity 
field which are in the form of potential vorticity. Thus, the bracket (^) may 
be interpreted as to define the Poisson bracket algebra of generalized potential 
vorticities on the flow space. We note finally that the bracket ([22] ) is induced 
from and is compatible with the Poisson bracket (|4]) on R x R 3 defined by 
the symplectic two-form Q. 
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4 Invariant differential forms 



We shall construct invariant differential forms of the velocity field associated 
with the infinitesimal symmetries. The techniques of constructing such in- 
variants together with the infinite dimensionality of the symmetry algebra 
will result in proliferation of invariant forms. However, a detailed analysis 
of them will lead us to identify a basic hierarchy to formulate the laws of 
general vorticity and helicity conservations. 

Since the time components £ k of symmetries ( ]T7| ) are conserved functions, 
the right hand side of the symmetry criterion (|i~3|) vanishes and this makes 
the three-form 

ot k = i(Uk) (/■*) = ikdx A dy A dz — (£ k v + Ufe) • dx A g?x A dt (25) 

which is independent of our choice for p^, an absolute invariant of dt + v. 
This means, 

C dt+V (a k ) = i([d t + v, U k ])(ji) = (26) 

where we used the identity 

C v i{U) - i(U)C v = i([V, U}) (27) 

and that v is divergence-free, a k s are closed via conservation of £ k s. The 
commutativity of the Lie derivative with the interior product for dt+v implies 
that the two-forms 

Ok = i(d t + v )(«fc) = — Ufc ■ <ix A dx — (u& x v) ■ dx A dt (28) 

are also absolutely invariant. They are closed by Eq.(^) and the closure of 
afc's or, by direct computation using Eqs.(|l^) and divergence-free property 
of Mfc's. Employing the Poincare lemma, we introduce one- forms 9k 

Q k = -d6 k , 9 k = ip k dt + A fc • rfx , u fc = V x A fc (29) 

and ipk are determined from the equations 

^ + u fc x v = W> fc . (30) 

Commutativity of the Lie derivative with the exterior derivative and the 
Poincare lemma imply that 9 k s are relative invariants. That means, their 
Lie derivatives with respect to dt + v are exact differentials 

C 9t+V (9 k ) = di(d t + v)(9 k ) = dxk , Xk = ipk + A fc ■ v (31) 

which follow from the definitions (^) and ( |29| ) . 
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4.1 A plethora of invariants 



It can be seen using the identities (13) and ( p?D that for an invariant form 
a and an infinitesimal symmetry Ui of v, the Lie derivative Cu t (a) and the 
interior product i(Ui)(a) are also invariants. In fact, is obtained from a k 
by taking interior product with dt + v which is trivially a symmetry. 

Similarly, we have the absolutely invariant two-forms a\ k = i(Ui)(ak) 
obtained from a^'s. Their Lie derivatives give daik = i([Ui, U k ])(fx) which 
may be regarded as redefinitions of a k s in the algebra of vector fields fll6|) . 

From 0^ = aok we get absolute invariants Oik = i(Ui)(Q k ) and dQik = 
i(dt + v)(dotik), the latter of which have similar interpretations of being re- 
definitions of Bfc's. 

The relatively invariant one-forms 9k result in conserved functions 9\ k = 
(i^ + Uj-Ajt of the velocity field whenever \k in Eqs.(|31|) are conserved under 
the flow of the symmetry U\. Lie derivatives of 9 k s add the differential —d9i k 
to the absolute invariants 0^ to produce new invariant one-forms which are 
also absolute invariants if 9\ k are conserved functions of v. 

In addition to taking Lie derivative and interior product with infinitesimal 
symmetries, forming wedge products of invariants also produces new invari- 
ants. For example, 6(, A 9j, 9 k A ai and 9 k A 0/ A 9i are relative invariants 
while & k A Bi = VJfe, I = 1, 2, 3, ... 

Thus, the above procedures seem to generate a plethora of invariant forms 
for the velocity field. However, we find some of them to be identical 

9 k A on = 9 lk fx = [Xkii + A fc • uj]/i (32) 
9 k AQ l A9 l = Xi(A fc ■ (33) 

and that they consist of functions available from more simple invariants. 
Moreover, we can get rid of the proliferation of invariants by imposing phys- 
ical restrictions and utilizing the Lie algebraic structure of the infinitesimal 
symmetries. 

We first observe that the invariants ak consist of two parts which decouple 
the velocity field from its symmetries and hence are physically not much 
relevant. To see this we write ak as the sum of two absolute invariants 

otk = £,kO - u fc • dx A dx A dt , a = T\ i i=l (dx l - v l dt) (34) 

where a is the well-known invariant three-form that appears in the formula- 
tion of continuity equation and has no effect in the definition of 0&'s. 
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The invariants a^, on the other hand, have a similar structure as those 
of involving appropriate differences of the symmetries u\ and Uk and thus 
can be taken to be equivalent. Likewise, due to the Lie bracket relations 
between symmetries we can regard da^ and dQ ik to be equivalent to (f25|) 
and (^|), respectively. 

5 General conservation laws 

This analysis of invariants leads us to conclude that appropriate expressions 
of the laws of general vorticity and helicity conservations may be best given 
in terms of the basic hierarchy consisting of G^ and 9k- We shall formulate 
the conservation laws as a restatement of infinitesimal conditions on these 
invariants in terms of integration over submanifolds of the flow domain. 

Proposition 3 Let S t and D t be two and three dimensional regions in RxR 3 
advected by the graph (t,g t ) of the flow of the velocity field v. Denote their 
smooth boundaries by dSt and dD t . Then, 

(1) The exact two-forms Qk = —d9k are the basic ingredients of general 
vorticity conservation law which can be expressed as 

tJ ®* = ~tA e k = ~i d Xk = 0, k = 0,1,2,... (35) 

(It J St dt J a St J dSt 

where Go is the symplectic two-form Q. 

(2) The relatively invariant three-forms 9kAd9i define the general helicity 
conservation laws 

^ / e k Ade t =f Xkd9i = -l dxkAO l = (36) 

at JD t JdD t JdD t 

if either ui is tangent to the projection onto spatial domain of dD t or, Xk = 
i{d t + v)(9k) = constant. 

(3) If, on the other hand, u\ is tangent to the surfaces Xk = constant 
then the generalized helicity densities in 9k A d9\ are Lagrangian invariants, 
that is, conserved functions of the velocity field. Moreover, they belong to the 
Poisson bracket algebra j\22] ) on R 3 of generalized potential vorticities. 

For an equivalent formulation, using the one- forms 9k, of the first part of 
proposition^) one replaces exactness and absolute invariance with Eqs.(pO|) 
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and relative invariance, respectively. In obtaining the second part, we used 
the properties of the velocity field that characterize it as the generator of 
volume preserving diffeomorphisms. Namely, it is divergence-free and is tan- 
gential to the boundary of D t . Eqs.(|35|) and the condition that the helicity 
is conserved for Xk — constant follow from the fact that the integral of exact 
forms over boundaries is identically (d o d ■ =) zero via Stokes' theorem. 
Thus, in the latter case, there will be no boundary terms contributed by U/, 
meaning that the three-forms are absolutely invariant and they result in the 
last conclusion of the proposition. To complete the proof of proposition(|3|) we 
shall supply, in the next section, coordinate expressions for these geometric 
arguments as well as for the global formalism of conservation laws. 

We end this section with a remark on the boundaries of the integration 
domains S t and D t which are two and three dimensional subsets of the time- 
extended space R x R 3 . For such decomposition of space-time, the one and 
two dimensional boundaries dS t and dD t do not close up, but rather have 
extention along a closed, finite interval of the time axis. Their projection 
onto spatial domain are a closed curve and a closed surface in the usual 
sense. For example, a circular helix whose finite-lenght axis is directed along 
the time axis is the boundary of a two-surface and thus is a closed curve in 
R x R 2 . 

6 Demonstrations 

We shall give explicit coordinate expressions for the general conservation laws 
of proposition(3). In the particular case of the Euler flow we shall show that 
the dynamical equations themselves, Bernoulli's equation, the Kelvin's cir- 
culation theorem, the Helmholtz' vorticity theorem, conservations of helicity 
and potential vorticity can be obtained from the invariance conditions of the 
basic hierarchy O^, 9^ and from the global formulations and (PS|). 

6.1 Bernoulli's equation 

For ideal fluids and for k = 0, we consider the invariance condition ( PT| ) of 
the canonical one-form in Eq.(||) with ip as given in Eqs.(^2|). Using the 
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rotational form 

<9v „,P 1 9 . .„„. 

__ VXW= _V(- + -^) (37) 

of the Euler equations we obtain from Eq.(pTI) the Bernoulli's equation 



|(^ )+v .V(^ + ^) = (38) 

where P + pv 2 /2 is the total (or stagnation) pressure. This says that the 
non-invariance of the kinetic energy under the flow of velocity field is due to 
the existence of pressure. 

Although, Eq.(^) is a direct consequence of (0), our purpose to in- 
clude this and the next example is to demonstrate the interplay between the 
present geometric framework for the Lagrangian description and the Eulerian 
dynamical equations. 



6.2 The Euler equations 

The hierarchy of Eqs.fl3"5|) contains, in particular, the dynamical equations 
for the Euler flow. This follows from the relative invariance of the one-forms 
9^. Eq.fl3"T|) for k = can be solved for the Lie derivative 

C dt+V (v ■ dx) = ^- t (\v 2 )dt + V(^ 2 ~^)-dx (39) 

where we used Bernoulli's equation to simplify the right hand side. Note that 
the Lie derivative is the one on R x R 3 . In terms of the differential operators 
on R 3 this expression reduces to 

^(v • dx) + C v {v • dx) = d(±v 2 - -) (40) 

which is the invariant form of the Euler equations of ideal fluid. It is this 
form that can be generalized for flows on general Riemannian manifolds of 
arbitrary dimensions 0. 



6.3 Helmholtz' vorticity theorems 

We first remark that in the integration of the two-forms 0^ over two-surfaces 
advected by the fluid motion, the term (u& x v) ■ dx A dt has no contribution 
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because its pull-back to Lagrangian coordinates by the solution g t of ([1|) van- 
ishes. Then, fl35|) becomes an expression for the general form of Helmholtz' 
vorticity theorem 

T + I ®k = -j: I u fe • dx A dx = (41) 
at J s t at J s t 

which reduces for k = 0, uo = w = V x u to the conservation of vorticity 
flux and, for uq = B to the Alfven's theorem of dynamo theory. 



6.4 Kelvin's circulation theorems 

For two-surfaces S t with smooth boundary dS t the invariance of the two- 
forms 0fc gives 

= ^/ d K = ll ( f A ^ rfx ( 42 ) 
at JdSt at JdSt 

via Stokes' theorem and the generic form of symmetry generators as curl 
vectors. This is a generalization of the Kelvin's circulation theorem for which 
the usual form is obtained when k = 0. It is the conservation of circulation 
of the velocity field Ao = v for the Euler flow and of the vector potential 
Aq = A for magnetization problems. 



6.5 General vorticity conservations 

An immediate generalization to the infinite dimensional left Lie algebra of 
symmetries of the well-known equivalence of the Helmholtz' and the Kelvin's 



theorems p[ follows directly from Eqs.(|35|), (|4"l|) and fl42|). For two-surfaces 



S t = gt(S) with smooth boundary 8S t we obtain 

/ ft*(ufc • dx A <ix) = / Ufc ■ <ix A <ix = f A^ • <ix = constant (43) 

JS JSt JdSt 

where gt is a flow of the velocity field. The constancy of the first integral is 
the preservation of coadjoint orbits of G^'s, of the second is the Helmholtz' 
theorem on conservation of flux of flow lines of the infinitesimal symmetry 
Ufc, and of the third is the Kelvin's theorem on conservation of circulation of 
vector potential A& for u^.. The fact that, for fixed u^, these expressions hold 
for all solutions gt of the velocity field is a manifestation of the right invariance 
on the cotangent bundle of the group of volume preserving diffeomorphisms. 
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6.6 General helicity conservations 

Eulerian invariants are conservation laws of dynamical equations (H) which 
are divergence expressions of the form 

dT 

— + V • P = (44) 

where the conserved density T and the flux P are functions of t, x, the Eu- 
lerian fields and their derivatives. Observing that Eq.([44"D can be expressed 
as the closure of the three-form 

- Tdx A dy A dz + P ■ dx A dx A dt (45) 

we shall derive the Eulerian conservation laws as secondary invariants ob- 
tained from the basic hierarchy and (P9"l). These are the relative in- 
variants 9k A dOi which are closed via Eqs.(pT8[) and (j30|). The corresponding 
conservation laws are expressed by Eqs. (|36"D of proposition(H). 

The case k = I = is exceptional because Q is non-degenerate. In this 
case, the identity 

d{9 Afi) + fiAfi = (46) 
gives the conservation of usual (magnetic) helicity 

|-(A • B) + V ■ ((A • B)v - (tjj + ip + v • A)B) = (47) 
at 

at which a hierarchy of others is anchored. For k, I ^ we find 

^-(A k ■ V x A,) + V • (A fc x (v x (V x A,)) - ^ fc V xA ( ) = (48) 
at 

associated with a pair of symmetries. Using the identity 

A fe x (v x u,) = (u, • A fc )v - (v • A fe )ui (49) 
and the divergence-free properties of v, we can write Eq.(£H) as 

d 

— (A k ■ u,) + v • V(A fe ■ u,) - u ; • V X k = (50) 
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from which one can immediately conclude the results in the second and the 
third items of proposition(3). Thus, the functions 

A • V x A , A • Vif x Vh , (pB-Vh , 
(uiXVxj^A-Vx^xVxj'-'A, k,l = 2,3,... (51) 

are the first few and the generic members of the conserved densities associated 
to the general helicity conservation (|36|) . The usual helicity and the potential 
vorticity are included as the first and the third elements of this hierarchy. 

The conservation of potential vorticity can also be obtained from the 
absolute invariance of the four-form Q A Q by integration over a region in 
R x R 3 . This, in turn, implies the constancy of the integral of the three- 
form 9 A Q over closed three dimensional domains in R x R 3 . For k = I = 
1,2,3,... the densities fl5"T|) vanish because the symmetries have Clebsch 
representations. 



6.7 Topological interpretations 

In the particular case of I = and k — 0, 1,2, the hierarchy uo ■ A& of 
invariants is related to the linking number of trajectories of vorticity uq = w 



or magnetic Uo = B fields with those of symmetries [pfl-||12|. For k = 0, 
the integral of this function coincides either with the self-linking number of 
vorticity or magnetic field, which is also called helicity, or, with the linking 
number of w and B expressed by the integrals of w • A or v • B. This 



latter result, when first proven with the second integrand in [12| for ideal 
magnetohydrodynamic equations, was interpreted to be unexpected because 
w is not frozen. Here, we see that it is sufficient to have B, and in general 
Ufc, to be frozen. 

The closed two-forms 0& can be associated with two-dimensional foli- 
ations of four-space by the integral curves of and V. The conditions 
0^ A 0; = means that intersection of leaves of any two such foliations is 
a one-dimensional foliation. The leaves of the latter are trajectories of the 
suspended velocity field. Moreover, for any k, I, (V,Uk,Ui) form a three- 
dimensional foliation. In this case, the integrals of the relative invariants 
9i A 9k A 0fc are shown to be related to the average linking number of the 
foliation of with the vector field Um defined by i(Uki)(fi) = d(9f- A 9{) 
1^| , ||14| . The characteristic form of this is the vector field = Xk^-i — Xi^k 



on R 3 
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We conclude from these examples that it is the hierarchy of secondary 
and higher order invariants of particle relabelling symmetry which contains 
topological informations about the flow domain. 

6.8 Lagrangian invariants 

Lagrangian invariants are functions which are conserved under the flow of 
the velocity field (|1|). The Eulerian conservation laws (0) reduce to the 
invariance of the functions u& • A;, k, I 7^ 0, 1 in the Lagrangian description 
when the condition u*. • Vxz = is satisfied. For the special case of I = for 
which xo — if > + f + v ' A, this can be verified using Eqs. fllOD and the identity 
v x (V x A) + (v • V)A = v l VAi. Similarly, one obtains from Eq.(|G]) and 
the symmetry condition for that ■ v is a Lagrangian invariant if 

u k -(F + ^Vv 2 )=0, k = 0,1,2,... (52) 

is satisfied by a general force field F. It can be checked that for the Euler 
equations of ideal fluid this coincides with the general criterion. 

The hierarchy of Lagrangian invariants includes the helicities v ■ w, A • B, 
the potential vorticities w • V/t, B ■ V/t as well as the invariants w ■ A and 
v • B of magnetization problems and, furthermore v ■ Vy? x V/t, v • Vy? x V/t 
which are obtained for the first two values of k, I. It can be seen using the 
Hamiltonian form with (|22|) of w^'s that the Lagrangian invariants • A; 
belong to the Poisson bracket algebra ( p2"D of generalized potential vorticities 
hk and can be obtained from the Poisson brackets of functions hk and hi. 

An immediate application of these conserved functions of the velocity 
field may be the construction of volume preserving diffeomorphisms which 
are also automorhisms of the symplectic two-form (|8|). This can be achieved 
by choosing a set y l (t, x) , i = 1, 2, 3 of functionally independent Lagrangian 
invariants which further satisfy the condition 

det = Vy 1 ■ Vy 2 x Vy 3 = 1 (53) 

for volume preservation. In particular, if we let y l = (p then Vy 2 x Vy 3 must 
be B (recall that we set p v = —B(lp) = 1). Thus, the other two Lagrangian 
invariants y 2 and y 3 are the Clebsch variables p, q for A (cf. Eqs.(|10])). From 
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proposition(l) we conclude that this particular diffeomorphism is the one 
which brings the symplectic two-form Q into canonical form. It follows from 
Eq.fl53|) and (f^) that 

{q,p} v = i- (54) 



7 Discussions and conclusions 

We presented a geometric approach to the explicit construction of symmetries 
starting from an initial one and utilizing a symplectic structure both of which 
are implicit in the Eulerian equations constraining the velocity field. We may 
conclude that, as the simplicity and the Lie-Poisson structure of Eulerian 
description derives from the particle relabelling symmetry 0-0, the Eulerian 
description, in its turn, enables us to construct the Lie algebra of these 
symmetries. 

Therefore, depending on the structure of the Eulerian dynamical equa- 
tions, the present construction may be generalized to other hydrodynamic 
systems. We demonstrated that the Euler flow in arbitrary Riemannian 
manifolds allows this sort of generalizations. A treatment for compressible 
flows may be based on the recognition that the so-called exponential vector 
field H exp(J V • vdt)w is a time-dependent infinitesimal symmetry of the 
velocity field. We refer to Refs. [{PI^ and [[14] for other examples of hydro- 



dynamic systems which admit a vectorial frozen-in field and hence can be 
analysed with the present geometric framework. 

It would be reasonable to compare the overall picture for invariants of 
three dimensional flow with that of two dimensional motion because, al- 
though, there is no essential difference in their group theoretical descrip- 
tions, they have been observed to have some different quantitative aspects 
[0-[|l6|,|l2|],|§,||. We showed that some properties of two dimensional flows 



such as the representation of area preserving diffeomorphisms as Hamiltonian 
vector fields and the existence of infinitely many enstrophy type integrals, 
can be realized, with appropriate modifications of geometric tools, for the 
three dimensional flows as well. 

Moreover, due to dimensional reasons we obtained a much richer structure 
of invariants for the present case. Namely, the geometric framework at our 
disposal led us to discriminate the invariants into three classes each having a 
certain degree of relation to the particle relabelling symmetry. Thus, for three 
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dimensional flows, in addition to general vorticity conservation, we can speak 
of a general helicity conservation and moreover, under certain conditions, of 
general Lagrangian invariants. 

The geometric character and the algebraic consequences of our construc- 
tion and classification of symmetries and invariants distinguish the present 
work from the approaches that employ the analysis of defining equations 
as the main tool || , |l7j , [18|] . Such constructions require that the algebraic 



structure of the solution set consisting of symmetries and invariants be sep- 
arately treated. Instead, we gave a recursive construction of symmetries via 
the Jacobi identity and, moreover, realized them as the Hamiltonian auto- 
morphisms of the symplectic two-form. This ensures that the infinitesimal 
symmetries thus obtained constitute a Lie algebra because, although not ev- 
ery infinite set of vector fields form a Lie algebra, Hamiltonian vector fields 
do so. 
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